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Abstract 



We determine the sphaleron transition rate using real time lattice simulations of the 
classical system. An improved definition of the lattice topological charge allows us to 
obtain a more reliable estimate of the transition rate. For an SU(2) Yang-Mills-Higgs 
system in the broken phase we find the transition rate to be strongly suppressed, and 
we have observed no sphaleron transitions in the range of coupling constants used. For 
a pure SU(2) Yang-Mills system in large volumes the rate behaves as k(q;„T)'', with 
K slightly decreasing as the lattice spacing is reduced. If the lattice size is reduced to 
about twice the magnetic screening length, the rate is suppressed by finite-size effects, 
and K is approximately proportional to the lattice spacing. Our rate measurements are 
supplemented by analysis of gauge field correlation functions in the Coulomb gauge. 
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1 Introduction 



Knowledge of the high-temperature baryon-number violation rate in the electroweak 
theory may be a key to origins of baryon matter. For temperatures close to and above 
the electroweak phase transition this essentially nonperturbative quantity can only be 
determined under two simplifying assumptions that (a) fermion degrees of freedom can 
be taken into account through modification of coupling constants, and that (h) the 
classical approximation is valid. The rate is then found by numerically following the 
real-time evolution of the Chern-Simons number in the classical (lattice) cutoff theory. 
This line of research was initiated in , and the last two years have seen much numerical 
effort to improve the simulations, as well as theoretical examination of the simplifying 
assumptions (mostly the second one). In our large-scale numerical simulation of the 
SU(2) Yang-Mills theory we found that, with the standard approximate discretization 
of topological charge density, the classical Chern-Simons number diffusion rate F is finite 
in the continuum limit: 

F = /^(a„T)^ (1) 

We also estimated the value of the dimensionless coefficient k in front of the law (P: 
K = 1.09 ± 0.05. Moore employed a different method for numerical determination of 



F, and his results agreed with ours. Tang and Smit |Q measured F in the SU(2) Yang- 
Mills-Higgs theory. In the low-temperature Higgs phase of the theory F is theoretically 
expected to be suppressed by the Boltzmann factor of the sphaleron potential barrier 
0, rendering Chern-Simons number diffusion inobservable in a realistic numerical sim- 
ulation. Surprisingly, no such strong suppression was found in Even more dramatic 
deviation from the sphaleron regime was found by Moore and Turok, raising suspicion 
that the residual diffusion in the Higgs phase is an artifact of the classical lattice the- 
ory 1^. Arnold, Son and Yaffe (ASY in the following), and later Arnold |p further 
questioned the validity of the classical result for F, arguing that the classical F must be 
proportional to the lattice spacing and therefore vanish in the continuum limit. More 
recently, Moore and Turok , employing an improved definition of Chern-Simons num- 



ber, found that the transition rate in the lower-temperature phase was zero within the 
measurement error, as expected from the sphaleron approximation. They also detected 
a dependence of the rate on the lattice spacing, but that dependence was considerably 
milder than the prediction p, |^ and left ample room for the existence of continuum 
limit. 

Our present work is an effort to further examine the classical result for F and to 
identify its dependence on the lattice discretization. Our new results show that, with the 
standard lattice definition of topological charge density, there are spurious contributions 
to F, directly related to the well-known difficulty with defining topological quantities in 
a compact lattice gauge theory. We introduce a new definition of the lattice topological 
charge, designed to strongly suppress these artifacts. With most of the artifacts removed, 
F is indeed inobservably low in the Higgs phase of the Yang-Mills-Higgs theory. In the 
pure Yang-Mills theory, the use of the standard definition of topological charge leads to 
a systematic error in F. Arnold, Son, and Yaffe conjectured that this systematic error 
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creates an illusion of a k independent of the lattice spacing, while the correct lattice rate 
vanishes in the continuum limit. We do find a tendency of the rate to decrease with the 
lattice spacing, but this tendency does not appear to be strong enough to rule out the 
existence of the continuum limit. There is nevertheless one instance in which the rate 
exhibits proportionality to the lattice spacing in a wide range of the latter: the Yang- 
Mills theory in a small volume, where finite-size effects give rise to a sphaleron barrier. 
In that case, however, there remains a possibility that the continuum rate exists, but is 
very small due to an exponential suppression of barrier configurations. 

The plan of this paper is as follows. We begin by explaining the cooling method for 
measuring lattice topological charge in Section |[ This method is used to place an upper 
bound on the rate in the low-temperature phase of the Yang-Mills-Higgs system (Section 
^ and to estimate the rate in the large-volume limit of the Yang-Mills theory (Section 
We also consider sphaleron transitions in the Yang- Mills theory in a small volume 
(Section A closer look at the ASY scenario is taken in Section H, where we study the 
dynamics of gauge fields in the Coulomb gauge. A discussion of the emerging picture is 
given in Section |^. 



2 Cooling method for the lattice topological charge 

Consider the standard lattice approximation for the topological charge per unit time 

^cs = ^ E (Eln + Ef_^,n) E Tr , (2) 



where 



Using the stndard notation, we assign to every lattice site j link matrices Uj^n in the 
fundamental representation of SU(2) and lattice analogs of color electric fields Ej^n, three 
of each, corresponding to the three positive lattice directions n. The usual plaquette 
variable is denoted by and a" are the three Pauli matrices. In fact, it is misleading 
to write the left-hand side of (@) as a time derivative of A^cs or any other fixed-time 
functional of the fields: the object on the right-hand side of (^ is not a total time 
derivative (TTD). This shortcoming of the lattice topological charge density (TCD) is 
well known and related to the difficulty in defining on a lattice an analog of the continuum 
A'cs with the correct properties under gauge transformations [|17|. Recently ASY M 



pointed out that deviations of the lattice topological charge per unit time from TTD 
may give rise to systematic errors in the transition rate. For smooth field configurations 
(0) tends to the continuum Nqs- However, thermal field configurations on the lattice are 
not smooth, hence (Q) does not approach a TTD in the continuum limit. 

The idea behind the cooling method is to replace the real-time trajectory in the phase 
space (obtained by integrating the classical equations of motion) by another path between 
the same endpoints, the one along which phase-space configurations are predominantly 
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smooth. Lattice topological charge, computed along such a path using the naive defini- 
tion @, would suffer much less from systematic errors than the one computed along the 
real-time trajectory. We shall first present the method for a wide class of Hamiltonian 
systems and later discuss its application to topological transitions. A simplified version 
of the method was used in [jrT], |T2| to study the time fiow of eigenvalues of the (lattice) 
Dirac Hamiltonian in the presence of the time dependent gauge field generated by the 
real time, finite temperature simulations, and it was shown that the cooling strips off 
the high energy modes, while it does not affect the (long wave length) eigenvalue modes 
which cross zero and which are responsible for the change of baryon number. The results 
reported in |1T[ provided the first proof that the sphaleron-like transitions observed in 



the real time lattice simulations are indeed long wave length excitations and not lattice 
artifacts. 

Let the real-time evolution of a dynamical system be described by a Hamiltonian 
depending on coordinates qi and momenta Pi {i = 1,...N): H = J^iPff^ + We 
introduce, along with the real time t, a cooling time r, of which all the dynamical 
variables are functions: Pi = Pi{t, t), = qi(t, r). The t dynamics at r = is given by 

dtq,{t, 0) = p,{t, 0); ^tp^{t, 0) = -d,y{qit, 0)), 

while the evolution in the r direction is an overdamped motion (cooling): 

dMt,r) = -d,yiqit,T)). (3) 

If we insist that the real-time equations for the coordinates dtqi{t,T) = piit.r) hold 
everywhere in the t, r plane, then the r dependence of momenta, consistent with (^, is 
given by 

d,p,{t,T) = -p,{T)dl^V{q{t,T)). (4) 

It is instructive to consider the r evolution described by (^-(^ in the vicinity of a static 
solution dq^V = 0. The equations of motion (§)-(§) can be linearized and dynamical 
variables can be chosen as eigenmodes of these linearized equations: drqi = —^^fqi, 
drPi = —^"iPi- Thus in the vicinity of a static solution cooling leads to an exponential 
decay of stable eigenmodes and exponential growth of unstable ones. Moreover, the rate 
of decay (growth) is especially rapid for high-frequency modes. 

Suppose there is a coordinate (call it Q) along which V is periodic. For simplicity 
we can assume that variations of Q are stable eigenmodes of the linearized equations of 
motion in the vicinity of periodic minima of V and the unstable ones in the vicinity of 
saddle points separating these minima. For sufficiently large r Q(t, 0) will be mapped by 
cooling to Q{t,T) in a close vicinity of a minimum, unless Q{t,0) lies close to a saddle 
point. As a result, if the Q{t, 0) motion consists of transitions between minima separated 
by fiuctuations close to a minimum, the picture will sharpen under cooling: fiuctuations 
of Q{t,T) near a minimum will have a smaller magnitude, whereas transitions between 
vacua (the time spent in the vicinity of a saddle point) will become shorter. If this 
indeed is the case for Q, two lessons can be learned. On the one hand, the average rate 
of transitions between minima for Q(t,T) should obviously be independent of r. On the 
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Figure 1: A H-shaped trajectory in the t, r plain used for determination of topological 
charge with cooling. The endpoints of a r = real-time trajectory, shown by a dotted 
line, are denoted by diamonds. 

other hand, care should be exercised if t is discretized for numerical purposes: for large 
r the transition times may become shorter than the t time step. 

Consider now cooling as a method for approximating topological charge. Lattice 
topological charge per unit time is not a TTD. For smooth field configurations a lattice 
TCD can be expanded in powers of the lattice spacing, the leading term in the expansion 
being the continuum TCD. For a thermal field configuration this expansion is formal since 
in three dimensions thermal configurations are not smooth. By cooling, however, these 
configurations can be made smooth enough for the expansion to make sense. Then a 
lattice TCD will give rise to much smaller deviations from a TTD. We are therefore led 
to the following method of determining the topological charge. For any two instances of 
real time consider evolution of lattice fields along a Il-shaped trajectory in the t, t plain, 
with r sufficiently large to give a smooth cooled configuration, while, at the same time, 
r -C t. Variation of topological charge with r is found by replacing t derivatives with r 
derivatives in the expression for TCD, whereas t derivatives at r 7^ 0, required for the 
horizontal part of the trajectory, are found using (^). If the space integral of a lattice 
TCD were a TTD, the topological charge along any such trajectory would be independent 
of r. This is not the case for an approximate lattice TCD. However, for sufficiently deep 
cooling the topological charge becomes nearly independent of r. Indeed, deviations from 
a TTD are suppressed along the horizontal part of the trajectory. The vertical parts of 
the trajectory, where deviations from a TTD are not suppressed, are of a fixed (short) 
duration and give rise to a small error in determination of topological charge. Unlike in 
the r = case, this error does not accumulate. We now write down the cooling equations 
for the Yang-Mills theory described on the lattice by the Kogut-Susskind Hamiltonian 



In doing so it should be kept in mind that link matrices U are SU(2) group elements, 
each belonging to its own group. Hence the cooling equation for a link U can be written 
using covariant derivatives D"" on the group to which U belongs: 




drU = -D'^HymD^'U. 



(5) 
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Here a is an adjoint SU(2) index. If we choose to work with right covariant derivatives, 
D'^U = —iUa°'. Since Hym is hnear in any given U, D'^Hym is obtained by retaining 
of Hym only terms proportional to U, with U replaced by D'^U. 

The cooling equations for the lattice electric fields E'^ are derived, in analogy with 
(^), by requiring, separately for each link, that the Hamiltonian equation of motion 
U = —iUa'^E'^ holds everywhere in the t, r plane. Denoting E = a°'E°' we then obtain 

drE = idriU^U) = {drU^)UE + iU^{HYM, drU}, (6) 

where {} are Poisson brackets and drU is substituted from (we refrain here from 
writing the somewhat cumbersome explicit expression for {HYM,drU}). 

The cooling method as described for the Yang-Mills theory applies essentially un- 
changed to the Yang-Mills-Higgs system whose Hamiltonian is 

Hh = Hym + E K.f + E \<l>^+n - uIM' + ^E (I'/'Z " ^'f ■ (7) 

j j,n j 

Our notation for the scalar fields is standard, i.e., to every lattice site j we assign 
the scalar doublet 0j along with its conjugate momentum Hj. For the gauge fields one 
simply replaces Hym by Hh in These equations are complemented by the cooling 

equations for the scalar fields and momenta, which are a straightforward implementation 
of ®-(®. 

3 The rate in the Yang-Mills-Higgs theory 

We used the cooling method in order to study the real-time evolution of Chern-Simons 
number in the SU(2) theory with and without the scalar field. In ([^) we chose A = 
0.09109, corresponding to the tree-level rriH/mw ratio of 0.423. Since the static proper- 
ties of the classical theory closely resemble those of the dimensionally reduced one, we 
expected our lattice system to undergo a strongly first order phase transition. This is 
indeed what we find: the system exhibits a metastability at beta = 11.25, with an easily 
detectable latent heat and a large discontinuity in the conventionally used "string bit" 
order parameter 

as shown in Figure ^ Therefore, according to the conventional wisdom |TB|, the sphaleron 
regime should set in at temperatures immediately below the transition. At /5 = 12 the 
sphaleron estimate for F involves the Boltzmann factor exp{—j3Esp^) which prohibits 
direct observation of sphaleron transitions in a feasible numerical simulation. 

Our numerical experiment at /5 = 12 reveals a dramatic difference between the naive 
and refined definitions of topological charge (Figure While the former is a random 
walk with deviations as large as 2 for a real-time trajectory of 2000 units, deviations 
of the latter never exceed 1. In fact, the motion of the refined Chern-Simons number 



6 



U 0.6 




1.14 



- 1.12 



11 

10 



1.08 



1.06 



Figure 2: The "string bit" order parameter (diamonds) and the average energy per degree 
of freedom normahzed by the energy of a free theory (pluses) for the Yang-Mills-Higgs 
theory with the tree- level m/^/miy = 0.423. Note the discontinuity of both quantities at 
(3 = 11.25. The error bars are smaller than the plotting symbols. The dotted lines are 
to guide the eye. 



is completely consistent with thermal fluctuations in the vicinity of one of the gauge- 
equivalent vacua and involves no diffusion between the vacua. In order to eliminate a 
possibility of diffusion with a time scale longer than 2000 units we followed the real-time 
evolution of the reflned A^^cs for 2 x 10^ time units. Again, we detected no sign of diffusion. 
Our results, together with observations of Ref . , strongly suggest that the slow diffusion 
of the naively deflned Chern-Simons is spurious. Using the refined definition instead, we 
can place an upper bound on the sphaleron transition rate, expressed in terms of the 
dimensionless coefficient k: k, = (7r/?)^/(W) < 0.01, where Vt is the total space-time 
volume of our simulation. Clearly, this upper bound cannot be used to confront the 
sphaleron estimate which it exceeds by many orders of magnitude. 



4 The rate in the Yang-Mills theory for small vol- 
umes 

Since we suspect the naive definition of topological charge density to give rise to sys- 
tematic errors in F, it is necessary to revise the results of for the Yang-Mills theory, 
obtained using that definition. While our ultimate goal is to determine F for lattice 



7 



A^cs -0-5 




2000 



Figure 3: Time history of Ncs in the Yang-Mills-Higgs theory with (the upper curve) 
and without (the lower curve) the coohng improvement. 



systems large enough to rule out finite-size effects, it is instructive to first consider 
smaller- volume systems. As we observed in P], finite-size effects in F set in as soon as 
the linear size L of the lattice is close to the inverse temperature /3. This phenomenon 
can be interpreted as follows. In the infinite-volume system there is no potential barrier 
separating topologically distinct vacua. Indeed, such a barrier is prohibited by the scale 
invariance of the classical theory. A finite size breaks the scale invariance and gives 
rise to the minimal barrier height of 18.15/L for the periodic boundary conditions used 
here [1^. Making use of this property, we can suitably adjust the jS/L ratio in order to 



study the A^cs behavior in the sphaleron regime. With this goal in mind, we performed 
a simulation at /3/L = 1. As Figure ^ shows, motion of the naively defined A'cs can 
be seen as consisting of three components: thermal fiuctuations, rapid transitions with 
a magnitude close to 1, and slow diffusion. If the refined definition is used, the slow 
diffusion component is absent from the A'cs time history. If we put this result together 
with our findings for the Yang-Mills-Higgs system, a consistent picture emerges: the 
true sphaleron transitions occur rapidly, while the slow diffusion of the naive A^cs is a 
discretization artifact. 

Next, we determine the transition rate F in the small- volume case. This is most easily 
and efficiently done by counting the clearly distinguishable sphaleron transitions. The 
rate is then the number of transitions per unit volume per unit time. The advantage 
of this procedure over determining the rate as the A^cs diffusion constant is that the 
sphaleron event count is the same with either A'cs measurement method. We therefore 
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Figure 4: Time history of iVcs in the Yang-Mills and without the cooling improvement 
at /3 = L = 12. Both curves exhibit rapid fluctuations and a sphaleron transition at 
t ~ 2000, but, for the naive definition of iVcs, there also is a slow diffusion upwards. 
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events 


111 
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73 


50 


time 


5.55 X 10^ 


8.2 X 10^ 1.23 X 10^ 


1.02 X 10^ 


1.44 X 10^ 


kappa 


0.195 ±0.018 


0.155 ±0.015 0.125 ±0.012 


0.112 ±0.013 


0.0676 ±0.0096 






Table 1: Effective k for (3 / L 


= 1. 




use the 


naive definition 


of A^cs in this part of the work. 







We have observed between 50 and 100 events ioi j3 = L = 10, 12, 14, 16, 20. Since 
any two consecutive events are widely separated in time and unlikely to be correlated, 
we determine the rate with the relative accuracy where n{(3) the total number 

of observed events for a given (3. Our findings, summarized in Table |I|, clearly show 
dependence of the effective k on (3, and hence on the lattice spacing. As can be seen in 
Figure |^, for 10 < /? < 16 this dependence is consistent with k oc but this simple 
proportionality cannot accommodate the /3 = 20 point in a convincing way. 

Clearly, our small-volume measurements fail to demonstrate the existence of contin- 
uum limit for the classical rate. Our simulations give support to ASY prediction k ~ 1//3 
except maybe for the largest value of (3. However, as we shall see in the following, a closer 
examination of the gauge field dynamics leaves room for alternative explanations of this 
discretization dependence of the small-volume rate. 
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Figure 5: Effective ft as a function of T = in the Yang- Mills theory for (3 = L. The 
dotted line is a linear fit through the origin using all the data except the /3 = 20 point. 



5 The large- volume rate in the Yang- Mills theory 

In order to eliminate the systematic error in F resulting from the naive definition of 
topological charge we performed a series of large-volume simulations similar to the one 
reported in but this time using the refined definition of A^^cs- The cooling method is 
very costly computationally. For this reason we were unable to produce a measurement 
sample of the same size as in with a comparable numerical effort. Our new results 
therefore contain larger statistical errors. In order to make use of the existing data from 
0] we adopted the following procedure. Whatever definition of Nqo, is used, the rate F 
is found by fitting the average squared topological charge {B^(t)) to a linear function of 
the time interval t in a suitable range of values of t. The rate is then the slope of that 
linear function. To rephrase, the rate can be determined given a numerical estimate for 
the joint probability distribution of Di = {B^{ti)), where (1 < z < /) is a discrete set 
of t values to be used in the fit. We assume this probability distribution to have the 
normal form 

Pi{D}) oc exp - A)(C-^).,(A - A)) , (8) 

where C is the estimated covariance matrix and is the estimate of obtained by 
averaging over the data sample {i.e., Di Di with increasing measurement accuracy 
|T3|). We shall take (|^) to represent results of our high-statistics study In the current 
simulation we measure the improved topological charge B'{t) along with the naive one 
B{t). We then find the corresponding Di and D-, which we can jointly call Aj, 1 < i < 21, 
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p 


T 

L 


length X number of trajectories 


K 


biased k 


10 


32 


2000 X 20 


1.48±0.08 


1.44±0.03 


12 


32 


2000 X 20 


1.36±0.09 


1.39±0.06 


14 


32 


2000 X 10 


1.33±0.12 


1.38±0.10 


14 


36 


1500 X 20 


1.27±0.15 




16 


36 


1500 X 20 


1.28±0.10 





Table 2: Summary of the large-volume study of k in the Yang-Mills theory. 

meaning that Aj = Dj for j < I and that Aj = D'^^i for j > I. We also estimate the 
covariance matrix C for {A}, leading to the joint probability distribution 11 for {A} 
completely analogous to (H): 

U{{D}, {D'}) cx exp (-1(A, - A,){C-%{A, - A,)) , (9) 

It is obvious that for a fixed set of {D} values of {D'j, averaged with the corresponding 
conditional probability, will be shifted away from {D } appearing in (P). This also holds 
true if the distribution (H) of {D} is given instead of a sharp set of values. We then 
average 11 over {D} with the weight P to obtain the corrected distribution for {D'}: 

P'{{D'}) = J [dD]U{{D}, {D'})P{{D}). (10) 

We refrain from writing down the somewhat cumbersome result of the trivial Gaussian 
integration ([101) . It is clear that P'{{D'}) is again a normal distribution, with corrections 
to both averages of D'- and to their covariance matrix. 

We now present results of our analysis. Table ^ summarizes the large-volume simu- 
lations we performed. Throughout the simulations the Hamiltonian time step was 0.05, 
and 5000 thermalization steps were performed before every Hamiltonian trajectory. 
We used cooling depth r = 0.2 for the improved determination of A^^cs- Within our 
measurement accuracy further cooling made no difference in the rate. As in our previous 
work, squared topological charge for a given time interval was first averaged within every 
Hamiltonian trajectory, and these averages, assumed independent, served as an input for 
further statistical analysis, as described. While the dimensionless coefficient k tends to 
decrease as a function of beta, our results cannot be said to favor the 1/(3 behavior of k, 
advocated in Ref. @]. 

6 Dynamics of gauge fields in the Coulomb gauge 

We have seen in both large and small volumes the tendency of the rate to decrease as 
the continuum limit is approached. A decrease of k as a function of (3 was advocated in 
Ref. Q. In a nutshell, Ref. argued as follows. Baryon- number violating processes 
mainly involve gauge fields with wavelengths ~ /3, whose motion is overdamped due to 



11 



interaction with hard modes. This overdamping leads to a behavior of k. While 
in the small-volume case we observed this behavior of k in a range of values of (3, the 
decrease in k for large volumes was significantly slower than ~ We therefore thought 
it useful to study the motion of the long-wavelength gauge fields in some detail. 

Gauge fields can only be defined by choice of a gauge. In this work we chose the 



Coulomb gauge: on the one hand, it is often used in the relevant analytical work ||T9 
on the other hand, lattice gauge-fixing techniques have been developed for this gauge. 
The lattice Coulomb gauge condition is usually defined for every lattice site j in terms 
of A,,„ = (l/(20(f/,,„-t/J, 



dAj = Y.iAj,n - A,^n,n) = 0, (11) 
n 

where n denotes the three lattice directions. For the SU(2) group this condition can be 
obtained by extremizing 

TrE^^.," (12) 

with respect to gauge transformations. This gauge is known to lead to Gribov copies; 
the unwanted copies are discarded by requiring the absolute maximum of (|12D. It is 
important to note that in a system with periodic boundary conditions (|TTp and (|T^ 
are preserved by global gauge transformations. This additional gauge freedom is not 
important as long as we only consider static properties of gauge fields, e.g., screening 
correlation functions, but it must be removed if we wish to compare gauge fields at 
different instances of dynamical evolution. We do so by imposing two additional global 
conditions: 

Trai Y: = Tra' E ^.M = Tra^ ^ A,,^ = 0. 
j j j 

In order to fix the gauge, we employed a combination of iterative overrelaxation and 
Fourier acceleration algorithms. The former was used to attain the accuracy of 10"'^ for 
Tii^dAy per site, and the latter to reduce Ti^dAy to below 10~^. We took Gribov copies 
into account by repeatedly performing first a random gauge transformation, then the 
Coulomb gauge fixing on the same field configuration, and choosing among the resulting 
gauge copies the one with the largest Tr J2j,n Uj^n- Within our measurement accuracy, 
this procedure had no effect on the observables. 

We first determined the magnetic screening length from screening correlation func- 
tions of the form 

G{z) = TiY.{A{z^)-A{z^ + z)), (13) 

where A{zq) is the average of A over the plane with one lattice coordinate fixed to Zq. 
We further average G{z) over the three directions of the displacement z (we will use the 
same notation for the average over directions). The magnetic mass is found by fitting 
G{z) to Cexp(-mL/2) cosh(m(z - L/2)). The results for (3 = 10, 12, 14 and L = 32, 
tabulated in Table |^, are accurately described by m = (2.630 ±0.005)//3 in lattice units. 
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/5 


10 


12 


14 


m 


0.27 ±0.02 


0.23 ±0.01 


0.184 ±0.006 


rriD 


0.712 ±0.025 


0.69 ±0.04 


0.68 ±0.04 



Table 3: Magnetic and Debye screening masses for the 32^ classical system. 



i.e., m = (0.657±0.001)5'^T in the continuum notation. Thus we find that in the classical 
theory and in the range of temperatures considered m increases with g'^T steeper than 
0.466g^T reported in Ref. for the high-temperature quantum theory. In that work 



correlation functions were computed at much lower P < 3.47 and in Landau, rather 
than Coulomb, gauge. Note, however, that both in Ref. [|I^ and in the current work 
G{z) could only be described by Cexp(— mL/2) cosh(m(2 — L/2)) for z close to L/2, 
whereas for small values of z the decay of G{z) is noticeably slower. That is to say, 
1/m underestimates distances for which transverse gauge fields are strongly correlated 
in space. 

Similarly, we estimated the Debye screening mass m/j from 

Gd{z) = Ti J2(E,{zo)E,{zo + z)), (14) 

where, e.g., E^^Zo) denotes the average of the Cartesian z component of the lattice color 
electric field over the lattice plane whose z coordinate is Zq. We further average Gj^ over 
the three lattice directions and fit the result to G exp{—moL/2) cosh(m£)(z — L/2)), just 
as we did for the magnetic mass. The results are summarized in the bottom line of Table 
^. Assuming m/5 oc l/\/^, our data is best described by = (2.328 ± 0.002) /^/J3, 
slightly above the perturbative estimate rriD = 2.01/v^ 0- Therefore for j3 < 14 
rriD/m < 3.7 in the classical theory, and these two scales are not completely separated. 
This point could be important for understanding our results for the rate. 

Next, we consider the dynamical behavior of gauge fields in the Coulomb gauge. 
We determined real-time correlation functions of fields with a definite latice momentum 
2'Kk/L in one of the three lattice directions, for < < 4. More specifically, we studied 

G{k,t)^Ke{Al^{t)Au,,{0))^ 

where p runs over the two transverse polarizations. We further averaged G{k, t) over 
the three directions of fc, and the notation G{k,t) is from now on used for that average. 
The resulting correlation functions, plotted in Figure exhibit the following important 
features. First the motion of the gauge fields is clearly overdamped for k = 0, 1: the 
corresponding G{k,t) decay without oscillations. For k = 3,4 the field motion is un- 
derdamped, and an oscillatory behavior of G{k,t) is clearly observed. The k = 2 case 
belongs to an intermediate regime, which, within our accuracy, we cannot classify as 
under damped or overdamped. 

Note that for k = 0,1 the G{k,t) /G{k,0) against t//?^ curves coincide (within our 
measurement accuracy) for all three values of (3. This coincidence, however, does not 
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hold for higher values of k. It appears that the higher the temperature, the higher is 
k at which the overdamping sets in. In the next section we confront these properties 
of the gauge field correlation functions with our rate measurements and with the ASY 
scenario. 

There is no obvious way to compare our G{0, t) with the W real-time correlation 
function computed in Ref. 0] for the Yang-Mills-Higgs theory. In that work the k = 
W decay rate was found to be very small and roughly proportional to 1/(3 in the 
high-temperature phase. However, that result was obtained using local objects which 
only approximate gauge fields (in the unitary gauge) sufficiently deep inside the low- 
temperature phase. 



7 Discussion 



We now summarize our findings. We saw that, with the spurious Nqs diffusion removed 
by the cooling method, the sphaleron rate in the low-temperature phase of the Yang- 
Mills-Higgs theory is too low to be measured in the space-time volume attained in our 
simulations. The best we can do in this case is place an upper bound on the rate. That 
upper bound exceeds the analytical sphaleron prediction by many orders of magnitude 
and hence it is of limited value. The rate in the low-temperature phase close to the 
phase transition is a crucial parameter in any scenario of electroweak baryogenesis. Our 



results, as well as those by Moore and Turok [10|, indicate that a straightforward real 



time simulation may not be an adequate tool for finding this, very possibly extremely 
small, rate. 

The focus of our numerical effort in this work is the approach of the classical rate 
to a continuum limit in the Yang-Mills theory. Our results for large volumes [L > 2/3) 
show that the dimensionless coefficient k tends to decrease with the lattice spacing. This 
decrease, however, does not appear to be rapid enough to match the ASY scenario, 
in which overdamping of gauge fields at the magnetic length scale leads to the zero 
continuum limit for n. In this respect our results are in a qualitative agreement with 



those of Moore and Turok |10], who used a different procedure for measuring topological 
charge. On the other hand, our values of k are systematically lower than those of Moore 
and Turok, except for the point closest to the continuum limit, f3 = 16. 

We observe a much stronger dependence of k on the lattice spacing in the small- 
volume setting P = L. Over a range of values of /?, this lattice spacing dependence is a 
simple proportionality. 

Our measurements leave open three possibilities for the behavior of the transition 
rate close to the continuum limit. 



1. The rate vanishes in the continuum limit, regardless of the lattice size to the 
magnetic screening length ratio. For small volumes this scenario is consistent with 
our observations. For large volumes, it can be reconciled with our data if a more 
rapid decrease in k sets in at larger /?. 
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2. Both the large- volume and the small- volume rates are finite in the continuum limit. 
This scenario is favored by our large-volume data. It could be that an asymptotic 
approach of the small-volume n to the continuum limit occurs at (3 values larger 
than those considered here. 

3. The large- volume rate remains finite, but the small- volume rate vanishes in the 
continuum limitQ. This possibility is favored by our data. 

Consider now the additional information on the gauge field dynamics obtained from 
the Coulomb-gauge correlation functions. As we have seen, motion of gauge fields with 
momenta and 27r/32 is overdamped. A closer examination shows, however, that this is 
not quite the overdamping of long-wavelength modes predicted by ASY. Indeed the time 
scale for the correlation function decay is proportional to for these wavelengths, if (3 
is varied while keeping the wavelength fixed. ASY predict that the decay time at a fixed 
wavelength should behave as 1/ (3. The decay times are short, of the order of several time 
units, i.e., shorter than the corresponding wavelengths, whereas, in the ASY derivation, 
it is essential that the decay times be much longer than the wavelengths. Moreover, as 
our magnetic mass measurement shows, the overdamping only occurs for modes with 
wavelengths somewhat above the magnetic screening length. Harder modes, presumably 
corresponding to the magnetic scale, are not overdamped. 

We can now discuss the three possibilities we have listed for the rate behavior in the 
continuum in the light of our knowledge of the gauge field dynamics in the Coulomb 
gauge. Consider sphaleron transitions in a small- volume Yang-Mills theory. With peri- 
odic boundary conditions, the sphaleron solution extends over the entire volume of the 



system Therefore it has a significant zero- momentum component of the gauge field 
whose motion is overdamped, with the decay time ~ (3"^ rather than (3 = L. It is 
possible that this overdamping gives rise to the 1/(3 behavior of k in small volumes. On 
the other hand, it appears unlikely that very soft overdamped modes play a role in the 
Nqo, diffusion in large volumes. Indeed, for our L = 32 lattices, these modes correspond 
to excitations whose linear size is close to L. If these modes were important for forma- 
tion of sphaleron- like configurations, we would expect a strong finite-size dependence of 
the rate. Such a dependence, however, is not observed (c/the /3 = 14 rate for L = 32 
and for L = 36). Hence the sphaleron-like transitions, whose spatial size is determined 
by the magnetic mass, should not be strongly affected by damping of these soft modes. 
Modes with momenta on the magnetic scale or harder are not overdamped. Thus it 
appears that the ASY mechanism does not operate in topology-changing processes in 
our large-volume simulations. 

The argument just given supports the third possibility on our list. That possibilty, 
however, is quite unusual: it requires that k vanishes above some critical value of (3/L. 
A way to avoid this situation could be as follows. Topology-changing transitions in 



^It seems very likely that k is a growing function of the volume. Thus we discard here yet another 
possibility, whereby the large- volume rate vanishes in the continuum, but the small- volume one remains 
finite. 
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small volumes could proceed through barrier configurations whose size is considerably 
smaller than L. These configurations need not have overdamped Fourier components. 
Indeed, all but the constant modes on our (3 = L lattices are not overdamped. These 
smaller barrier configurations are stronger suppressed by Boltzmann factors than the 
sphaleron, but that does not necessarily preclude their finite contribution to the rate in 
the continuum. That continuum, /^-independent k might be so small that, in the range 
of f3 values considered, k is still dominated by the l/jS term. This is a way in which the 
second possibility on our list may be realized. 

Finally, the first possibility on our list cannot be ruled out. While it is not favored 
by our data, the case may be otherwise in a weaker-coupling regime. In the current work 
the Debye to magnetic mass ratio is below 3.7, so these two scales may not be completely 
separated. It may be necessary to further increase f3 in order to reach the range of validity 
of perturbation theory on which the ASY scenario is substantially based. Since m/mu 
only decreases as 1/a/5, while the rate decreases at least as rapidly as exploring 
couplings much weaker than those presently considered is likely to be computationally 
costly. 
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Figure 6: Real-time correlation functions G{k^ t) of gauge fields in the Coulomb gauge 
plotted against t/ P^. The value of k is indicated on each plot. Each G{k, t) is normalized 
by its respective value at t — 0. Diamonds correspond to (3 = 10, plus signs to j3 — 12 
and squares to P — lA for a 32^ lattice. In addition, the open circles in the k — plot 
show ^(0, t)/G{Q, 0) for /? = 14 on a 14^ lattice. Wherever the errorbars are not shown, 
they are smaller than the plotting symbols. 
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